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GEOMETRY OF MEMBRANES 

Z. C. TU 

Communicated by Boris Konopeltchenko 



Abstract. This review reports some theoretical results on the Geometry of mem- 
branes. The governing equations to describe equilibrium configurations of lipid 
vesicles, lipid membranes with free edges, and chiral lipid membranes are derived 
from the variation of free energies of these structures. Some analytic solutions to 
these equations and their corresponding configurations are also shown. 



1. Introduction 

Membranes are very crucial to living organisms. They are the barriers of cells 
and ensure cells to be relatively isolated individuals but still able to exchange 
some materials between the inner sides and outer surroundings through specific 
ways due to the fancy properties of membranes. Membranes usually consist of 
lipid bilayers mosaicked various kinds of proteins. They are also the key factors 
to determine shapes of some kinds cells. In particular, the biconcave discoidal 
shape of red blood cells is regarded as a result of minimizing the free energy of 
membranes under the area and volume constraints E1QTI because red blood cells 
have no complex inner structures. The equilibrium configurations of membranes 
have attracted much attention of mathematicians and physicists [ 19.32, 34, 3 811461 . 
A membrane is thought of as a 2-dimensional (2D) smooth surface in Euclidean 
space E 3 because its thickness is much smaller than its lateral dimension. The 
first step to investigate configurations of membranes is constructing a free en- 
ergy functional by consideration of symmetry. Then the governing equations to 
describe the equilibrium configurations can be derived by variation of the free en- 
ergy with some constraints. The next task is seeking for solutions to satisfy the 
governing equations and comparing the results with typical experiments. 

In this review, we will present some purely theoretical results on Geometry of 
membranes. For simplicity, we merely focus on structures of lipid membranes 
and only select the theoretical problems that both physicists and mathematicians 
are interested in. The governing equations to describe equilibrium configurations 
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of lipid structures are derived. Several solutions to these equations and their cor- 
responding geometries are also investigated. The rest of this review is organized 
as follows. In Sec. |2l we give a brief introduction to preliminary in mathemat- 
ics and physics including surface theory and variational method based on moving 
frame, and Helfrich's model of lipid bilayer. In Sec. [3] we derive a shape equation 
that describes equilibrium configurations of lipid vesicles — closed lipid bilay- 
ers. Then we discuss some analytic solutions and their corresponding configura- 
tions which include surfaces of constant mean curvature, torus, biconcave discoid, 
cylinder-like vesicles, and so on. In Sec.lU we investigate a lipid membrane with 
free edge(s). The shape equation and boundary conditions describing equilibrium 
configurations of the membrane are derived. Then we discuss the compatibility 
between the shape equation and boundary conditions, and verify five theorems of 
non-existence. In Sec. [5j we construct the free energy functional of chiral lipid 
membranes in terms of symmetric argument and then derive the governing equa- 
tions to describe their equilibrium configurations by variational method. Some 
analytic solutions and their corresponding configurations are also shown. In the 
last section, we give a brief summary and a list of related open questions. 

2. Preliminary in Mathematics and Physics 

In order to continue our discussion, we first introduce several key mathematical 
and physical concepts and tools that will be used in the following sections. 

2.1. Surface Theory Based on Moving Frame 
2.1.1. Moving Frame Method 

A membrane can be regarded as a smooth orientable surface embedded in E 3 . The 
properties of surface such as mean curvature and Gaussian curvature determine 
the shape of membrane. As shown in Fig. [Q each point on surface M. can be 
represented a position vector r. At that point we construct three unit orthonormal 
vectors ei, and e3 with e3 being the normal vector of surface M. at point r. 
The set of right-handed orthonormal triple- vectors {ei, 62,63} is called a frame 
at point r. Different points on the surface have different vectors r, ei, e2, and e3, 
thus the set {r; ei, e2, 63} is called a moving frame. 

Let us imagine a mass point that moves from position r to its neighbor position r' 
on the surface. The length of the path is denoted by As. Then we can define the 
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Figure 1. The moving frame on a surface. 



differentiation of the frame as 



and 



dr = lim (r' - r) = o;iei + w 2 e2, (1) 
As->-0 



dei = lim (e- - e*) = Wj,-e,-, (i = 1, 2, 3) (2) 

As— >0 



where lo\, lu2, and coij = 1, 2, 3) are 1-forms, and 'd' is the exterior differ- 
ential operator [7]. The 1-form u>ij is ant-symmetric with respect to i and j, that 
is ujij = — ujji. Here and in the following contents without special statements, the 
repeated subscripts represent summation from 1 to 3. Additionally, the structure 
equations of the surface can be expressed as Q : 

dwi = U)\2 A U)2 

dui 2 = ^21 A uji (3) 
duij = u ik A u kj = 1,2,3), 

and 



^23 / V " C / V W 2 

where 'A' represents the wedge production between two differential forms. The 
matrix ( ° ^ j is the representation matrix of the curvature tensor. Its trace 
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and determinant are two invariants under the coordinate rotation around e 3 which 
are denoted by 

2H = a + c and K = ac - b 2 . (5) 

H and K are called the mean curvature and Gaussian curvature, respectively. 
They determine the shape of the surface and can be expressed as 2H = — (1/Ri + 
I/R2) and K = I/R1R2 with two principal curvature radii R\ and R2 at each 
point on the surface. 

Now consider a curve on surface M.. Its tangent vector is denoted by t. Let (f) be 
the angle between t and ei at the same point. Then the geodesic curvature k g , the 
geodesic torsion r g , and the normal curvature k n along the direction of t can be 
expressed as B2l : 

kg = (d<f> + wi 2 )/ds 

T g = b cos 2cj) + (c — a) cos 4> sin 4> (6) 
k n = a cos 2 cj) + 2b cos 4> sin <f> + c sin 2 </>, 

where ds is the arc length element along t. If t aligns with ei, then <j) = 0, 
kg = wi2/ds, T g = b, and k n = a. In the principal frame, the geodesic torsion 
and normal curvature can be expressed as 

cos 2 4> sm2 4> 

K = 5 5 — , T g = (1/Ri - I/R2) cos (psm (p. (7) 

Ki -ft 2 



2.1.2. Stokes' theorem and related identities 



Stokes' theorem is a crucial theorem in differential geometry, which reads 



to = / dui, (8) 

where D is a domain with boundary d®. w is a differential form on dD. In 
particular, f^dto = for a closed domain 23. 

From Stokes' theorem, we can derive several identities listed as follows 114211461 . 
(i) For smooth functions / and h on 2D subdomain D C M. , 

/ (fd*dh - hd*df) = I (f*dh - h*df) 
[ (fd*dh - hd*df) = I (f*dh - h*df) (9) 
[ (fdidh - hdidf) = I (f*dh - h*df), 
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where * is Hodge star operator satisfying *u>i = lo 2 and = d and 

* are generalized differential operator and generalized Hodge star which satisfy 
df = /1W13 + /2^23 and id/ = f x u 2 z - f 2(^13 if df = fxU\ + / 2 ^2 B21|43|. 
(ii) If u and v are two vector fields defined on 2D subdomain T> C .M, then [44J 



where the 'dot' represents the inner product of vectors. For simplicity, Eqs.® 
and (fTOl) are still called Stokes' theorem in this review. They are widely used in 
the variational process. 

Additionally, we can also define the gradient, curl, divergent, Laplace operators, 
etc. on the surface in terms of the differential operators and Hodge stars. They are 
summarized as follows M451I46II : 



(V xu)dA = d(u • dr), (V ■u)dA = d(*u • dr) 
(V -u)dA = d(iu • dr), (V ■u)dA = d(*u • dr) 
V/ • dr = df, V/ • dr = df 
(V 2 /) dA = d* df, (V • V/) dA = d * df 
(V • V/) dA = did/ 

(u-V/)d^ = u-drA*d/ (11) 

(u • V/) dA = u • dr A *df 

(u • V/) dA = u • dr A id/ 

(V 2 u) dA = d * du, (Vu) • dr = du 

(Vu: Vv)dA = duAdv 

S ■ dr = — CJ12, (V • S)dA = —d * W12, 



where dA = wi A W2 and S are the area element and spin connection of the sur- 
face, respectively. A represents calculating dot production and wedge production 
simultaneously. 




(u • *dv — v • *du) 




(10) 
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2.2. Helfrich's Model 

In 1973, Helfrich proposed a spontaneous curvature model to describe the free 
energy of lipid membranes by analogy with a bent box of liquid crystal in Smectic 
A phase iPTO . The free energy is in fact a functional defined in the space of shapes 
of membranes, which reads 



F H 



M 



h {2 H + co) 2 + kK 



dA, (12) 



where M. , H , and K represent the membrane surface, mean curvature, and Gaus- 
sian curvature, respectively. k c and k are two bending moduli. The former should 
be positive, while the latter can be negative or positive for lipid membranes. Co 
is called spontaneous curvature, which reflects the asymmetrically chemical or 
physical factors between two leaves of lipid bilayers. 

On the other hand, the spontaneous curvature model can also be obtained from 
symmetric argument. A lipid membrane can be locally regarded as 2D isotropic 
elastic entity. Thus the local free energy density / should be invariant under ro- 
tational transformation around the normal direction of the membrane surface. In 
other words, it should be a function of H and K because H and K are the funda- 
mental invariants of the surface under rotational transformation. Up to the second 
order terms of curvatures, it can be expanded as 

f c = A + A X H + A 2 H 2 + A 3 K, (13) 

which can be rewritten as 

f c = ^(2H + c f + ~kK (14) 

with omitting an unimportant constant. This is no other than the integrand in 
Eq. (fT2l) . Therefore, the spontaneous curvature model is of general significance 
not only for lipid membranes, but also for other membranes consisting of isotropic 
materials. The following discussions are mainly based on Helfrich's spontaneous 
curvature model. 



2.3. Variational Method Based on Moving Frame 

To obtain governing equations that describe equilibrium configurations of lipid 
membranes, we need to minimize the free energy. That is, we should calculate 
the variation of functionals defined in the space of shapes of membranes. In this 
space, each shape might be expressed as a function, but the definitional domain 
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is not a fixed planar domain. It is difficult to deal with this case by using the 
traditional calculus of variation. Now we will introduce the variational method 
based on the moving frame developed by the present author and Ou-Yang 1141 11421 
l46l . which is available to deal with variational problems on surfaces. 

Any infinitesimal deformation of a surface can be achieved by a displacement 
vector 

5r = v = $7,ej (15) 

at each point on the surface, where 5 can be understood as a variational operator. 
The frame is also changed because of the deformation of the surface, which is 
denoted as 

5ei = Ojj-ej (i = 1,2,3), (16) 

where Qij = —Qji = 1, 2, 3) corresponds to the rotation of the frame due to 
the deformation of the surface. From 5dr = ddr, ddej = ddej, and Eqs. CD)-©, 
we can derive 

5u>i = dv • ei - u>2^2i = dtti + fl 2 ^2i + ^3^31 - w 2 J7 2 i 
5u)2 = dv • e 2 - u>\tti2 = dfl 2 + ^1^12 + ^3^32 - coi^u (17) 
StUij = dflij + QuUJij — LOuVtij 
dv • e 3 = d0 3 + VL\uj Vi + $7 2 u; 2 3 = + Sl 2 3^ 2 - 

These equations are the essential equations of the variational method based on the 
moving frame. With them as well as Eqs. (0]) and ([5]), we can easily derive 

5dA = (divv - 2Hn 3 )dA 
8{2H) = [V 2 + (AH 2 - 2K)]U 3 + V{2H) • v (18) 
SK = V • V0 3 + 2KHn 3 + VK-v. 

Using the above equations Q-©, (|9l)- (|TTb . (fTTT ) and (fT8T >. we can deal with almost 
all variational problems on surfaces. 

3. Lipid Vesicles 

Most of lipid molecules are amphiphiles with a hydrophilic head group and two 
hydrophobic hydrocarbon tails. When a quantity of lipid molecules disperse in 
water, they will assemble themselves into a bilayer vesicle as depicted in Fig. [2] 
due to hydrophobic forces. In this section, we will theoretically understand vari- 
ous configurations of lipid vesicles. 
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Figure 2. A cartoon of lipid vesicle. 



3.1. Shape Equation to Describe Equilibrium Configurations 

A lipid vesicle can be represented as a closed surface Ai. Its equilibrium shape 
is determined by minimizing Helfrich's free energy (fT2)) under the constraints of 
constant area and volume because experiments reveal that the area of lipid mem- 
branes are almost incompressible and the membranes are impermeable for the 
solutions in both sides of the membranes. Thus we can introduce two Lagrange 
multipliers A and p to replace these two constraints, and then minimize the fol- 
lowing functional 



where M., A and V represent the membrane surface, total area of the vesicle and 
volume enclosed by the vesicle. A and p can also be understood as the apparent 
surface tension and osmotic pressure of the lipid vesicle. 

The Euler-Lagrange equation of functional (fT9l ) can be derived by using the vari- 
ational method presented in Sec. |2j which reads 1127 1128 [l42l 



with reduced parameters p = p/k c and A = X/k c . This formula is called the 
shape equation because it describes the equilibrium shapes of lipid vesicles and 
represents the force balance along the normal direction of membrane surfaces. 

Now let us consider an axisymmetric vesicle generated by a planar curve shown 
in Fig. [3] Rotate the curve around z— axis and then mirror with respect to the hor- 




p - 2XH + (2H + c )(2H 2 - c H - 2K) + 2V 2 H = 



(20) 
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Figure 3. The generation curve for an axisymmetric vesicle. 



izontal plane, p is the revolving radius, and ip is the angle between the tangent of 
the curve and the horizontal plane. That is, dz/dp = tan ip. Through simple cal- 
culations, we can obtain —2H = h = s'mi/j/p + (sin-^)', K = simjj(simjj)' / p, 
and V 2 {2H) = — (cos V'/p) (p cos iph')' , where the 'prime' represents the deriva- 
tive with respect to p. Substituting these relations into Eq. (1201 ). we derive 



The equivalent form of this equation is first derived by Hu and Ou-Yang lfT4l . The 
shape equation (f2TT > of axisymmetric vesicles is a third-order differential equation. 
Following Zheng and Liu's work 1531 . we can transform it into a second order 
differential equation 



with an integral constant 770- It is found that the shape equation of axisymmetric 
vesicles obtained by Hu and Ou-Yang degenerates into that derived by Seifert et 
al. 11371 when 770 = in Eq. (1221 or for vesicles with spherical topology free of 
singular points 1133 11551 . 

3.2. Analytic Solutions and Corresponding Configurations 

Now we will show several analytic solutions to shape equations (1201 or (1221 of 
lipid vesicles that we have known till now. They correspond to surfaces of constant 
mean curvature, torus, biconcave discoid, and so on. Most of these solutions are 
found by Ou-Yang and his coworkers. We will see that only sphere, torus, and 
biconcave discoid can correspond to lipid vesicles. We also recommend the reader 




(21) 



cos iph' + (h — cq) sin t/jif)' — A tan ip + 



2?7q — pp 2 tan if) 



(h - c ) 2 = (22) 



2pcos^ 2 
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to note the classic paper by Konopelchenko [17 ] who expressed Eq. [20] in terms 
of the inverse mean curvature and density of squared mean curvature, and then 
found a broad variety of its solutions and corresponding possible shapes. 

3.2.1. Surfaces of Constant Mean Curvature 

Obviously, surfaces of constant mean curvature (including sphere, cylinder, and 
unduloid shown in Fig. [4]) can satisfy the shape equation. First, let us consider a 
spherical surface with radius R. Then H = — 1/R and K = 1/R 2 . Substituting 
them into Eq. (1201 . we derive 

pR 2 + 2XR- c (2 - c R) = 0. (23) 

This equation gives the relation between sphere radius R, spontaneous curvature 
cq, reduced osmotic pressure p, and reduced surface tension A. 




Figure 4. Surfaces of constant mean curvature: sphere (left), cylinder (mid- 
dle), and unduloid (right). 

Next, H = —1/2R and K = for a cylindrical surface with radius R. Then 
shape equation d20l requires 

2pR 3 + 2XR 2 - 1 + c 2 R 2 = 0. (24) 

For other surfaces of constant mean curvature such as unduloid, H is a constant 
but K is not a constant. In terms of shape equation ([20l . we obtain 

H = -c /2 (25) 

and p = —Acq. 
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Note that cylindrical surface and unduloid are in fact not closed surfaces. Alexan- 
drov proved that "An embedded surface (no self-intersection) with constant mean 
curvature in E 3 must be a spherical surface" [1]. Thus we can only observe one 
kind of vesicles of constant mean curvature — sphere. 

3.2.2. Torus 

As shown in Fig. [5j A torus is a revolution surface generated by a circle with 
radius r rotating around an axis in the same plane of the circle. The revolving 
radius R should be larger than r. The torus can be expressed as vector form 
{(R + r cos <p) cos 9,(R + r cos p) sin 9, r sin p>}. Through simple calculations, 
we have 2H = — (R + 2r cos p)/r(R + r cos ip) and K = cos pjr(R + r cos </?). 
Substituting them into Eq. (|20T >. we derive 

[(2cgr 2 - 4c r + 4Ar 2 + 2pr 3 )/u 3 ] cos 3 ip 
+ [(5cor 2 - 8c r + lOAr 2 + 6pr 3 )/i/ 2 ] cos 2 tp 
+ [(4cor 2 - 4c r + 8Ar 2 + 6pr 3 )/u] cos tp 
+2/v 2 + {c 2 r 2 - 1) + 2(pr + A)r 2 = (26) 

with v = R/r. If v is finite, then Eq. (1261 ) holds if and only if the coefficients of 
{1, cos p, cos 2 cp, cos 3 p} vanish. It follows 2Ar = co(4 — c$r), pr 2 = — 2co and 

v = Rjr = y/2. (27) 

That is, there exists a lipid torus with the ratio of its two generation radii being 
y/2 (called V2 torus by Ou-Yang 11291 ). which was confirmed in the experiment 
lf2ll . It is also found that nonaxisymmetric tori |[36l constructed from conformal 
transformations of v2 torus also satisfy the shape equation. 




Figure 5. Torus (left) and its generation curve (right). 



To check the consistency, we need also verify that the above toroidal solution can 
indeed satisfy shape equation (1221 of axisymmetric vesicles. It is not hard to see 
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y/2 torus can be generated from a curve expressed as 



sin i/> = (p/r) ± s/2. 



(28) 



Substituting it into Eq. (|22l . we arrive at 2Xr = co(4 — cqt), pr 2 = —2cq and 
7/o = -1/r ^ 0. 

3.2.3. Biconcave Discoid 

For < cqPb < e, the parameter equation 



corresponds to a planar curve shown in Fig. [6] Substituting it into Eq. (l22l) . we 
have p = 0, A = 0, and rj = 2cq / 0. That is, a biconcave discoid generated 
by revolving this planar curve around z-axis can satisfy the shape equation of 
vesicles. This result can give a good explanation to the shape of human red blood 
cells under normal physiological conditions Il9l l22ll24l . 



A small comment is that rjo = 2co ^ reflects the singularity at two poles of 
the biconcave discoid. Whether does this singularity exist in real red blood cells? 
What is the biological meaning of this singularity? Or does there exist a normal 
solution to the shape equation that can also explain the shape of red blood cells? 
These open questions need further discussions. 

3.2.4. Unduloid-Like and Cylinder-Like Surfaces 




(29) 




Figure 6. Biconcave discoid (left) and its generation curve (right). 



If we only concern solutions to the shape equation, two cases are also widely 
discussed. The first case M201I23I is an axisymmetric surface generated by planar 
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curve satisfying 

sinV , = _L_ (_£_ + PHL\ _ * _ 2) (0 < p mCo < 4/3). (30) 

The generated surface abides by the shape equation with p = ~2copf n , A = 
2/ ' p 2 m — Cq/2, and rjo = 2cq — "ijcQp 2 m . This surface has the unduloid-like shape 
but with nonconstant mean curvature. 

The second one is a cylinder-like surface generated by a planar curve translating 
along the normal of the plane. If we denote the curvature of the curve as k, then 
the geometric quantities of the generated surface can be expressed as 2H = — k, 
K = 0, and V 2 {2H) = —k, where the 'dot' above k represents the derivative 
with respect to the arc length of the curve. Thus shape equation (|20l degenerates 
into |2l|l9j|54l: 

p + \K + K 3 /2-k = (31) 
with A = A + Cq/2. The above equation is integrable, which results in 

k 2 = £,o + 2pK + Xk 2 - k 4 /4 (32) 

with an integral constant £q- This equation can be further solved in terms of 
Elliptic functions Il2l l49ll54ll56l . 

It is necessary to note that these two cases do not correspond to real vesicles 
because they are not closed surfaces. 

4. Lipid Membranes with Free Edges 

The opening-up process of liposomal membranes by talin ll35l was observed, 
which gives rise to the study of equilibrium equation and boundary conditions 
of lipid membranes with free exposed edges. This problem was theoretically in- 
vestigated by Capovilla et al. [5 ] and Tu et al. iHTI in terms of different methods. 
In this section, we will present these theoretical results and subsequent advance- 
ments. 

4.1. Shape Equation and Boundary Conditions to Describe Equilibrium Con- 
figurations 

As shown in Fig. |7J a lipid membrane with a free edge can be expressed as an 
open smooth surface {M) with a boundary curve (C) in geometry. Because the 
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Figure 7. An open smooth surface {M) with a boundary curve (C). 
{ei,e2,e3} forms the frame at some point on the surface. {t,b,e3} also 
forms the right-handed frame for the point in C such that t is the tangent of 
C and b points to the side that the surface is located in. 



free exposed edge is energetically unfavorable, we assign the line tension (energy 
cost per length) to be 7 > 0. Then the free energy functional that we need to 
minimize can be expressed as 



F 



M 



-(2H + coY + kK 



dA + XA + jL, 



(33) 



where L is the total length of the free edge. 

By using the variational method introduced in Sec. |2j we can arrive at a shape 
equation [41] 

(2H + c )(2H 2 - c H - 2K) - 2XH + V 2 {2H) 

and three boundary conditions [41 J 



0, 



(34) 



\{2H + c ) + kn T 



c 



0. 



[-2dH/db + 7^n + ~kt g 



c 



0, 



[(1/2) (2H + c ) 2 + kK + A + 7^] 



c 



0. 



(35) 
(36) 
(37) 



where A = X/k c , k = k/k c , and 7 = j/k c are the reduced surface tension, re- 
duced bending modulus, and reduced line tension, respectively. K n , K g , and r g 
are the normal curvature, geodesic curvature, and geodesic torsion of the bound- 
ary curve, respectively. The 'dot' represents the derivative with respect to the arc 
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length of the edge. Equation d34l expresses the normal force balance of the mem- 
brane. Equations (l35T)-(l37l) represent the force and moment balances at each point 
in curve C. Thus, in general, the above four equations are independent of each 
other and available for an open membrane with several edges. 



Figure 8. Outline of an axisymmetric open surface. Each open surface can 
be generated by a planar curve AC rotating around z axis, i/j is the angle 
between the tangent line and the horizontal plane. 

Now we consider axisymmetric membranes. When a planar curve AC shown 
in FigJH revolves around z axis, an axisymmetric surface is generated. Let tp 
represent the angle between the tangent line and the horizontal plane. Each point 
in the surface can be expressed as vector r = {pcos <ft, psin 0, z(p)}, where p 
and <f) are the radius and azimuth angle that the point corresponds to. Introduce a 
notation a such that a = 1 if t is parallel to dr / d(j), and cr = — 1 if t is antiparallel 
to dr/dcj) in the boundary curve generated by point C. The above equations (l34l)- 
(l37l ) are transformed into 



z 



A 




C 



P 




(38) 
















(39) 



(40) 



(41) 
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with h = smi/j/p + (simp)' and K = sin^(sin?/>)'//9. The 'prime' represents 
the derivative with respect to p. 

Shape equation (I38T ) is integrable, which reduces to a second order differential 
equation 

cos ibti + (h-co) sin iM' - A tan ib + V ° , - (h - c ) 2 = (42) 

p cos ip 2 

with an integral constant rjo B71 . This equation is equivalent to Eq. (1221) with 
zero osmotic pressure. The configuration of an axisymmetric open lipid mem- 
brane should satisfy shape equation (l42l and boundary conditions d39T>— (l4Tb . In 
particular, the points in the boundary curve should satisfy not only the boundary 
conditions, but also shape equation (l42l because they also locate in the surface. 
That is, Eqs. (l39l)- (|4TI) and ((42]) should be compatible with each other in the edge. 
Substituting Eqs. d39T)-(l4~l"l) into (|42~1 ). we derive the compatibility condition WT\ 
to be 

Vo = 0. (43) 

It is a necessary (not sufficient) condition for existence of axisymmetric open 
membranes. Under this condition, the shape equation is reduced to 

cosiph' + (h - c ) sin^' - A tan tp ^- (h - c ) 2 = 0, (44) 

while three boundary conditions are reduced to two equations, i.e. Eqs. (l39l > and 
(SD- 



4.2. Theorems of Non-Existence 

Now our task is to find analytic solutions that satisfy both the shape equation and 
boundary conditions. An obvious but trivial one is a circular disk with radius R. 
In this case, Eqs. (l34l-(l37T) degenerate to 

\R + fy = 0. (45) 

Can we find nontrivial analytic solutions? We will prove several theorems of 
non-existence in this subsection, which imply that it is almost hopeless to find 
nontrivial analytic solutions. 

Theorem 1 There is no open membrane being a part of a spherical vesicle. 
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Proof. For a sphere with radius R, we can calculate H = —1/R, n n = —1/R 
and T g = in terms of Eq. © because a = c = — 1/R and b = for a sphere. 
Boundary condition (|36l ) cannot be abided by. Thus an open membrane cannot be 
a part of a spherical vesicle. 

Theorem 2 There is no open membrane being a part of a cylindrical surface. 

Proof. For any line element on the surface of a cylinder with radius R, we can 
calculate K n = — cos 2 9/R from Eq. © where 9 is the angle between the line 
element and the circumferential direction. Additionally, H = — is a constant. 
If k = 0, then boundary condition (l36l ) results in K n = 0, that is 9 = tt/2. The 
line along this direction is not a closed curve, and so cannot be as an edge of a 
membrane. If k ^ 0, then boundary condition (I35T ) results in n n = (cq — l/R)/k, 
which implies 9 should be a constant. The unique closed curve is a circle, i.e. 
9 = and n n = —1/R. But r g = if 9 = 0, then contradicts with boundary 
condition (l36l) . Thus an open membrane cannot be a part of a spherical surface. 

Theorem 3 There is no open membrane being a part of a curved surface with 
constant mean curvature. 

Proof. Two special surfaces (sphere and cylinder) with constant mean curvature 
are discussed in the above theorems. Now we only need to investigate surfaces 
with constant H but nonconstant K. From shape equation of open membranes, 
we derive two possible cases: (i) H = —cq/2 ^ and A = 0; (ii) H = cq = 
and A / 0. 

In the former case, if k = 0, then boundary conditions (l36l > and (l37l) result in 
k„ = Kg = 0. Thus the curvature of boundary curve is k = y / 1 n 2 n + k 2 , = 0. That 

is, this curve is a straight line which is not closed curve. If k ^ 0, then boundary 
conditions (|35T > and (l36l) results in n n = and r 9 = constant. Using Eq. (Q, 
we derive the principal curvatures for the points in the curve are constant. Then 
Eq. (l37l) requires n g = constant. That is, the curvature and torsion are constant 
in the curve. The unique closed curve is a circle. But r g = for a circle. Let c\ 
and C2 represent the two principal curvatures, 9 is the angle between the tangent 
of the curve and one principal direction at each point in the curve. Then we have 
two equations: K n = c\ cos 2 9 + C2 sin 2 9 = and r g = (c2 — c\) sin 9 cos 9 = 0. 
Substituting these two equations and K = c\Ci into Eq. (1371 ). we obtain K g = 0. 

Then k = ^Jk^ + k 2 = 0, which contradicts with the preassumption of a circle. 
In the latter case, H = cq = 0, similar to the proof in the former one, it also 
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leads to a contradiction. Thus there is no open membrane being a part of a curved 
surface with constant mean curvature. 

Theorem 4 There is no open membrane being a part of Willmore surface. 

Proof. Let us consider the scaling transformation r — > Ar, where the vector r 
represents the position of each point in the membrane and A is a scaling parameter 
0. Under this transformation, we have A — > A 2 A, L — > AL, H — > A~ l H, and 
K — > A~ 2 K. Thus, Eq. (1331) is transformed into 

F(A) = [ \{k c /2){2H) 2 + kK]dA 

Jm 

+ 2k c c A [ HdA + (X + k c c 2 /2)A 2 A + jAL. (46) 
Jm 

The equilibrium configuration should satisfy dF/dA = when A = 1. Thus we 
obtain 

2c / HdA + (2A + cl)A + 7 L = 0. (47) 
Jm 

This equation is an additional constraint for open membranes. 

Willmore surfaces satisfy the special form of Eq. (1201) with vanishing A and Co 
I5T1 . Because 7L > 0, thus the constraint d47l) cannot be satisfied when A = 
and cq = 0. That is, there is no open membrane being a part of Willmore surface. 



Figure 9. Schematics of several impossible open membranes with free 
edges: parts of sphere, cylinder, unduloid, torus, biconcave discodal surface 
(from left to right). 

Corollary There is no open membrane being a part of v2 torus. 

Proof. Substituting v2 torus into Eq. (I34l . we obtain cq = and A = 0. That is, 
\[2 torus is a Willmore surface. In terms of theorem 01 we arrive at this corollary. 

Theorem 5 There is no axisymmetric open membrane being a part of a biconcave 
discodal surface generated by a planar curve expressed by sin i/j = —cqp hi(p/ Pb) 
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Proof. A biconcave discodal surface Il22ll24l generated by a planar curve ex- 
pressed by sin^ = — copln(p/ ps) with non-vanishing constants cq and ps- To 
avoid the singularity at two poles, we may dig two holes around the poles. Sub- 
stituting this equation into shape equation (l42l . we obtain A = and 770 = 2co- 
That is, the biconcave discodal surface can be a solution to the shape equation. 
However, 770 = 2co 7^ contradicts to compatibility condition (|43T ). Thus there is 
no axisymmetric open membrane being a part of this biconcave discodal surface. 

In Fig. |9j we show several impossible open membranes with free edges in terms 
of the above theorems. These theorems suggest that it is hopeless to find exactly 
analytic solutions to the shape equation and boundary conditions of open lipid 
membranes. Quasi-exact solutions or numerical simulations are highly appreci- 
ated. 

4.3. Quasi-Exact Solutions 

Here the quasi-exact solution is defined as a surface with free edge(s) such that 
the points on that surface exactly satisfy the shape equation, and most of points 
in the edge(s) abide by boundary conditions. In fact, the proves to theorems |2] 
and [3] implies two possible solutions as shown in Fig. [TOj One is a straight stripe 
along the axial direction of cylinder, another is a twist ribbon which is a part of a 
minimal surface (H=0). 




Figure 10. Schematics of two quasi-exact solutions: straight snipe along 
axial direction of cylinder (left) and twist ribbon which is a part of a minimal 
surface (right). 

Let we consider a long enough straight stripe along the axial direction of cylinder 
that satisfies shape equation (l34l ). that is, A = (1 — CqR 2 )/2R 2 . The long enough 
configuration ensure us to omit the boundary of two ends. The lateral edges are 
straight lines which have n n = n g = r g = 0. Thus boundary conditions (1331 ) and 
d36l ) are trivial. The third boundary condition results in A = (1 — cqR) 2 /2R 2 . 
Thus we arrive at A = and R = 1/cq. That is, a long enough straight stripe 
along the axial direction of cylinder with R = 1/co is a quasi-exact solution. 

Next, a twisted ribbon with pitch T and width 2uq can be expressed as vector form 

{u cos 99, u sin (p, atp} with \u\ < uq, \tp\ < 00 and |or| = T/27T. From simple 
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calculations, we have 

H = 0, K = -a 2 /(u 2 + a 2 ) 2 (48) 
for the points on the surface, and 

K n = 0,K = -a 2 /(u 2 + a 2 ) 2 

Kg = uo/(«o + « 2 ) (49) 
T g = a/(ul + a 2 ) 

for the points on the edges. 

It is easy to see that equation (l48l) can satisfy shape equation (l34l when Co = 0. 
Then equation d49l) naturally validates boundary conditions (1351 and (l36l) . The 
last boundary condition (1371 ) leads to A = [ka 2 — 71*0(^0 + a 2 )]/(«Q + a 2 ) 2 , 
which can be satisfied by proper parameters A, k, 7, uq and a. That is, the twist 
ribbon is indeed a quasi-exact solution. 

5. Chiral Lipid Membranes 

In fact, our above discussions only concern lipid membranes where lipid molecules 
are in Smectic A phase. In this phase, lipid molecules almost point to the normal 
direction of the membrane surface. However, there are also many kinds of chi- 
ral lipids in cell membranes. At body temperature, chiral lipids usually form 
Smectic C* phase. They are tilting from the normal direction in a constant angle. 
It is necessary to develop Helfrich's spontaneous curvature model introduced in 
Sec. 12.21 to cover the Smectic C* phase. Based on symmetric argument or Frank 
energy in the theory of liquid crystal, many theoretical models and results were 
achieved !!121ll6[|2^[30l[3T1,[39ll40l . These theoretical models contain much com- 
plicated terms and many parameters, which make it is impossible to derive the ex- 
act governing equations for describing equilibrium configurations of chiral lipid 
membranes. Here we will discuss a simplified version proposed by the present 
author and Seifert [45]. It is found that this concise theory can still explain most 
of experimental phenomena. 

5.1. Constructing the Free Energy 

The free energy density for a chiral lipid membrane are supposed to consist of the 
following contributions. 
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Figure 11. Right-handed orthonormal frame {ei, e2, 63} at any point in a 
surface where e3 is the normal vector of the surface, (a) Surface without 
boundary curve, (b) Surface with boundary curve where t is the tangent 
vector of the boundary curve, and b, in the tangent plane of the surface, is 
perpendicular to t. (Reprint from Ref. [45 1) 

(i) The bending energy per area is still taken as Helfrich's form (fT4l . That is, we 
neglect the anisotropic effect of lipid molecules' tilting on the bending moduli. 

(ii) The energy per area originating from the chirality of tilting molecules has the 
form 03 

fch = -hr m , (50) 

where h reflects the strength of molecular chirality. Without losing the generality, 
here we only discuss the case of h > 0. r m is the geodesic torsion along the unit 
vector m at each point. Here m represents the projected direction of the lipid 
molecules on the membrane surface. If we take a right-handed orthonormal frame 
{ei, e2, 63} as shown in Fig.QT] rn can be expressed as m = cos <pe\ + sin <pe2, 
where <p is the angle between m and ei. then geodesic torsion r m and normal 
curvature K m along m can be expressed as the similar form of Eq. (©. 

(iii) The energy per area due to the orientational variation is taken as 

/ (W = (A; / /2)[(Vxm) 2 + (V-m) 2 ] ; (51) 

where kf is a constant in the dimension of energy. This is the simplest term of 
energy cost due to tilting order invariant under the coordinate rotation around the 
normal of the membrane surface. By defining a spin connection field S such that 
V x S = K, one can derive (V x m) 2 + (V • m) 2 = (V0 — S) 2 through simple 
calculations 11251 . 

The total free energy density adopted in the present paper, G = f c + f c h + fov, 
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has the following concise form: 

G = ^-(2H + c ) 2 + ~kK-hT m + ^-v 2 , (52) 

with v = V</> — S. This special form might arguably be the most natural and con- 
cise construction including the bending, chirality and tilting order, for the given 
vector field m and normal vector field e3. 

5.2. Governing Equations to Describe Equilibrium Configurations 

The free energy for a closed chiral lipid vesicle may be expressed as 

F= f GdA + XA + pV, (53) 

Jm 

where A is the area of the membrane and V the volume enclosed by the vesicle. 
A and p are two multipliers to implement area and volume constraints. 

Using the variational method mentioned in Secf2j we can obtain two governing 
equations to describe equilibrium configurations B31 as 

2h(K m — H) — k f V 2 <f) = (54) 

and 

2X7 2 H + (2H + c ){2H 2 - c H - 2K) - 2XH + p 
+h[V • (mV x m) + V x (mV • m)] 

+k f [(K v - H)v 2 - Vv: Ve 3 ] =0 (55) 

with reduced parameters h = h/k c , kf = kf/k c , p = p/k c , and A = X/k c . K m 
and k v are the normal curvature along the directions of m and v, respectively. 
When writing Eq. (1541 . we have selected the proper gauge such that V • S = 0, or 
else V 2 (p should be replaced with V 2 </> — V • S. Additionally, we do not consider 
singular points for closed vesicles different from toroidal topology. 

Consider a chiral lipid membrane with a free edge as shown in Fig. [TTb . Its free 
energy can be expressed as 

F= f GdA + XA + jL, (56) 
Jm 

where A is the area of the membrane and L the total length of the edge. 7 repre- 
sents the line tension of the edge. 
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Using the variational method mentioned in Secf2l we can obtain the governing 
equations to describe equilibrium configurations of membrane surfaces as 

2h(K m — H) — k f V 2 (f) = (57) 

and 

2V 2 H + (2H + c )(2H 2 - c H - 2K) - 2\H 
+h[V ■ (mV x m) + V x (mV • m)] 

+k f [(«v - H) v 2 - Vv : Ve 3 ] = 0. (58) 

Simultaneously, the boundary conditions obeyed by the free edge are derived as 
US: 

v b = 0, (59) 
(l/2)(2# + c ) 2 + kK- hr m + (fc//2)v 2 + A + j Kg = 0, (60) 
(2H + c ) + kK n - (h/2) sin 20 = 0, (61) 
jK n + kt g — 2dH/db — h(v t + <j>) sin 2<p + kfK n vt = 0, (62) 

where n n , r g and K g are the normal curvature, geodesic torsion, and geodesic 
curvature of the boundary curve (i.e., the edge), respectively. v b and vt are the 
components of v in the directions of b and t, respectively. The 'dot' represents the 
derivative with respect to arc length parameter s. <ft is the angle between m and t at 
the boundary curve. Equations (T59l-(l62l) describe the force and moment balance 
relations in the edge. Thus they are also available for a chiral lipid membrane with 
several edges. 

5.3. Solutions and Corresponding Configurations 

Now we will present some analytic solutions to the governing equations of chiral 
lipid membranes. 

5.3.1. Sphere 

For spherical vesicles of chiral lipid molecules with radius R, r m is always van- 
ishing because a = c = 1/ R and 6 = 0. Thus the free energy d53l is independent 
of the molecular chirality and permits the same existence probability of left- and 
right-handed spherical vesicles. This is uninteresting case in practice. 
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5.3.2. Cylinder 

Here we consider a long enough cylinder with radius R such that its two ends can 
be neglected. The cylinder can be parameterized by two variables s and z which 
are the arc length along the circumferential direction and coordinate along axial 
direction, respectively. Let be the angle between m and the circumferential 
direction. Then Eqs. (1541 and (|55T > are transformed into B3I : 

k f (cl> ss + ^ zz ) + (h/R) cos 20 = 0, (63) 

and 

h[2(4> 2 z + (j> sz ) sin 2(f) + (<f> ss - </) zz + 4<^0 S ) cos 2cj>} + \/R 

+p + (c 2 - l/R 2 )/2R + k f [{<g - 4> 2 S )/2R + </> sz /R] = 0. (64) 

where the subscripts s and z represent the partial derivatives respect to s and z, 
respectively. 

It is not hard to see that 4> = 7r/4 and 2pR 3 + 2XR 2 - 1 + CqR 2 = can satisfy 
the above two equations. Thus a cylinder shown in Fig. [T2k with uniform tilting 
state (tilting angle (p = n/ 4) is a solution. 



(a) (b) 




Figure 12. Two possible chiral lipid membranes: (a) Cylinder with uniform 
tilting state; (b) Torus with uniform tilting state. 

5.3.3. Torus 

A torus is a revolution surface generated by a cycle with radius r rotating around 
an axis in the same plane of the cycle as shown in Fig. [5] It can be expressed as 
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vector form { (R + r cos tp) cos 8,(R + r cos tp) sin 9, r sin tp}. Equation (1541 ) is 
transformed into |j45ll : 



d 2 ct) d 
+ 



v + cos tp dO 2 d(p 



(z , + cosy ,)_ 



uhr 

cos 2(j) = 0, (65) 



A; 



where <fi is the angle between m and the latitude of the torus, while v = R/r is 
the ratio between two generated radii of the torus. 

The uniform tilting state (4> = — vr/4) satisfies Eq. (l65l) and makes — J hr m dA to 
take the minimum. With <p = — 7r/4, Eq. (|55T > is transformed into RBI : 



(2 - k f )/u 2 + (c^jr 2 - 1) + 2{pr + X)r z 

>T,„ i Q\„2 i (?- 3\ /,, 





— 4cor 




- 8c r 




— 4cor 



4^.r + lOAr 2 + 3k f + 6pr 3 )/V 2 ] cos 2 
4c r - 2hr + 4Ar 2 + 2fc/ + 2pr 3 )/u 3 } cos 3 = 0. (66) 

Because v is finite for a torus, then the above equation holds if and only if the 
coefficients of {1, cos tp, cos 2 tp, cos 3 tp} vanish. It follows that 2Ar 2 = (4rco — 
r 2 Cg) — 3k f + 2hr, pr 3 = 2k f — 2rco — hr and 



u = ^(2-k f )/(l-kf). (67) 

Thus a torus with uniform tilting state as shown in Fig. [T2b is an exact solution 
to governing equations of chiral lipid vesicles. The ratio of two generation radii 
satisfies Eq. (|67"T ). which increases with kf. Especially, v = y/2 for kf = 0, which 
leads to the v2 torus of non-tilting lipid molecules ||29l . Since this kind of torus 
was observed in the experiment [21 J, tori with v > v2 for < kf < 1 might also 
be observed in some experiments on chiral lipid membranes. 

5.3.4. Twisted Ribbons 

Here we consider a quasi-exact solution for the governing equations to describing 
equilibrium configurations of chiral lipid membranes with free edges. Two long 
enough twisted ribbons with lipid molecules in different tilting states are shown 
in Fig. [13] Similar to Sec 14.31 a twist ribbon can be expressed as vector form 
{u cos tp, u sin tp, aip} with \u\ < W/2, \tp\ < oo and |q| = T/2tt. Equation (I57T) 
is transformed into P3I : 



, , ucj) u + (j)^ 2hasm2(/) 
kf 4> uu H n - ; — 2 H 2~, — 2~ = °- ^ 68) 
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Figure 13. Long enough twisted ribbons with lipid molecules in different 
tilting states: (a) m is perpendicular to the edges; (b) m parallels the edges. 
Arrows represent the projected directions {m} of the tilting molecules on the 
ribbons' surface. (Reprint from Ref. J45|) 



where <f> is the angle between m and the horizontal. 

If we only consider the uniform tilting state, the above equation requires <f> = or 
7r/2. It is easy to see that = minimizes —hj T m dA for a < while = tt/2 
minimizes —hf T m cL4 for a > because T m = —a cos 2cj)/(u 2 + a 2 ) El. 
Thus we should take (j) = for a < and = tt/2 for a > 0. The former 
case corresponds to Fig. [T3h where m is perpendicular to the edges; the latter 
corresponds to Fig.[T3b where m is parallel to the edges. Both for 4> = and tt/2, 
Eq. (l58l) leads to cq = for non-vanishing a. Among the boundary conditions 
d59l>-(l62l). only Eq. (l60l) is nontrivial, which reduces to 



A(l + x 2 )a 2 -(h- yx)\a\ + ^ ■ „. = (69) 

2(1 + x z ) 

with x = W/2\a\. Solutions to this equation exists for proper parameters. Thus, 
there are indeed twist ribbons in two states as shown in Fig.[[3l they have different 
chirality and tilting angles. 
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6. Summary and Conclusion 

In the above discussions, we have presented some theoretical results on the Ge- 
ometry of membranes, which include the surface theory and variational method 
based on moving frame, the governing equations to describe equilibrium con- 
figurations of various lipid structures derived from the variation of free energy 
functionals, some analytic solutions to these equations and their corresponding 
configurations. We only focus on the pure theoretical researches and miss all ex- 
perimental and numerical results related to our topic on which gentle readers may 
consult Refs. ESSHMISISISSQ 

Although many theoretical advancements have been achieved, there are still a 
lot of challenges waiting for further investigations. Several key open questions 
among them are listed as follows. 

(i) Lipid vesicles of multi-components. Cell membranes contains many kinds of 
lipids. At body temperature, different kinds of lipids usually separate into several 
lipid domains. Lipid vesicles with two or several domains have been investigated 
from experimental and theoretical levels Oil 1011 151142114611501 ,53 1. However, there is 
still lack of strictly exact solutions to the governing equations B42N461 describing 
the vesicles with multi-domains. 

(ii) Other solutions on the shape equations of lipid membranes. We have only 
found a few analytic solutions to the governing equations of lipid structures. 
Whether are there other solutions, in particular to the simplest equations (1201 and 
(|22l ? Or can we prove that there is no other analytic form except the solutions 
that we have mentioned? 

(iii) Generalized boundary conditions for open lipid membranes. Although we 
have investigated the boundary conditions of lipid membranes with free edges, 
there are still other kind cases, such as confined edges, contact lines, and so on. 
Can we develop a generalized variational principle covering such cases? 

(iv) Non-orientable membranes. All membranes that we have considered are ori- 
entable membranes. How can we deal with the non-orientable membranes, such 
as Mobius band B71 ? 
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